Genetic underpinnings of arthropod community distributions in Populus trichocarpa.

Community genetics seeks to understand the mechanisms by which natural genetic variation in heritable host phenotypes can encompass assemblages of organisms such as bacteria, fungi, and many animals including arthropods. Prior studies that focused on plant genotypes have been unable to identify genes controlling community composition, a necessary step to predict ecosystem structure and function as underlying genes shift within plant populations. We surveyed arthropods within an association population of Populus trichocarpa in three common gardens to discover plant genes that contributed to arthropod community composition. We analyzed our surveys with traditional single-trait genome-wide association analysis (GWAS), multitrait GWAS, and functional networks built from a diverse set of plant phenotypes. Plant genotype was influential in structuring arthropod community composition among several garden sites. Candidate genes important for higher level organization of arthropod communities had broadly applicable functions, such as terpenoid biosynthesis and production of dsRNA binding proteins and protein kinases, which may be capable of targeting multiple arthropod species. We have demonstrated the ability to detect, in an uncontrolled environment, individual genes that are associated with the community assemblage of arthropods on a host plant, further enhancing our understanding of genetic mechanisms that impact ecosystem structure.

Validation of a metabolite-GWAS network for Populus trichocarpa family 1 UDP-glycosyltransferases.

Metabolite genome-wide association studies (mGWASs) are increasingly used to discover the genetic basis of target phenotypes in plants such as Populus trichocarpa, a biofuel feedstock and model woody plant species. Despite their growing importance in plant genetics and metabolomics, few mGWASs are experimentally validated. Here, we present a functional genomics workflow for validating mGWAS-predicted enzyme-substrate relationships. We focus on uridine diphosphate-glycosyltransferases (UGTs), a large family of enzymes that catalyze sugar transfer to a variety of plant secondary metabolites involved in defense, signaling, and lignification. Glycosylation influences physiological roles, localization within cells and tissues, and metabolic fates of these metabolites. UGTs have substantially expanded in P. trichocarpa, presenting a challenge for large-scale characterization. Using a high-throughput assay, we produced substrate acceptance profiles for 40 previously uncharacterized candidate enzymes. Assays confirmed 10 of 13 leaf mGWAS associations, and a focused metabolite screen demonstrated varying levels of substrate specificity among UGTs. A substrate binding model case study of UGT-23 rationalized observed enzyme activities and mGWAS associations, including glycosylation of trichocarpinene to produce trichocarpin, a major higher-order salicylate in P. trichocarpa. We identified UGTs putatively involved in lignan, flavonoid, salicylate, and phytohormone metabolism, with potential implications for cell wall biosynthesis, nitrogen uptake, and biotic and abiotic stress response that determine sustainable biomass crop production. Our results provide new support for in silico analyses and evidence-based guidance for in vivo functional characterization.

A k-mer based approach for classifying viruses without taxonomy identifies viral associations in human autism and plant microbiomes.

Viruses are an underrepresented taxa in the study and identification of microbiome constituents; however, they play an essential role in health, microbiome regulation, and transfer of genetic material. Only a few thousand viruses have been isolated, sequenced, and assigned a taxonomy, which limits the ability to identify and quantify viruses in the microbiome. Additionally, the vast diversity of viruses represents a challenge for classification, not only in constructing a viral taxonomy, but also in identifying similarities between a virus' genotype and its phenotype. However, the diversity of viral sequences can be leveraged to classify their sequences in metagenomic and metatranscriptomic samples, even if they do not have a taxonomy. To identify and quantify viruses in transcriptomic and genomic samples, we developed a dynamic programming algorithm for creating a classification tree out of 715,672 metagenome viruses. To create the classification tree, we clustered proportional similarity scores generated from the k-mer profiles of each of the metagenome viruses to create a database of metagenomic viruses. The resulting Kraken2 database of the metagenomic viruses can be found here: https://www.osti.gov/biblio/1615774 and is compatible with Kraken2. We then integrated the viral classification database with databases created with genomes from NCBI for use with ParaKraken (a parallelized version of Kraken provided in Supplemental Zip 1), a metagenomic/transcriptomic classifier. To illustrate the breadth of our utility for classifying metagenome viruses, we analyzed data from a plant metagenome study identifying genotypic and compartment specific differences between two Populus genotypes in three different compartments. We also identified a significant increase in abundance of eight viral sequences in post mortem brains in a human metatranscriptome study comparing Autism Spectrum Disorder patients and controls. We also show the potential accuracy for classifying viruses by utilizing both the JGI and NCBI viral databases to identify the uniqueness of viral sequences. Finally, we validate the accuracy of viral classification with NCBI databases containing viruses with taxonomy to identify pathogenic viruses in known COVID-19 and cassava brown streak virus infection samples. Our method represents the compulsory first step in better understanding the role of viruses in the microbiome by allowing for a more complete identification of sequences without taxonomy. Better classification of viruses will improve identifying associations between viruses and their hosts as well as viruses and other microbiome members. Despite the lack of taxonomy, this database of metagenomic viruses can be used with any tool that utilizes a taxonomy, such as Kraken, for accurate classification of viruses.

Genome-Wide Association Study of Wood Anatomical and Morphological Traits in Populus trichocarpa.

To understand the genetic mechanisms underlying wood anatomical and morphological traits in Populus trichocarpa, we used 869 unrelated genotypes from a common garden in Clatskanie, Oregon that were previously collected from across the distribution range in western North America. Using GEMMA mixed model analysis, we tested for the association of 25 phenotypic traits and nine multitrait combinations with 6.741 million SNPs covering the entire genome. Broad-sense trait heritabilities ranged from 0.117 to 0.477. Most traits were significantly correlated with geoclimatic variables suggesting a role of climate and geography in shaping the variation of this species. Fifty-seven SNPs from single trait GWAS and 11 SNPs from multitrait GWAS passed an FDR threshold of 0.05, leading to the identification of eight and seven nearby candidate genes, respectively. The percentage of phenotypic variance explained (PVE) by the significant SNPs for both single and multitrait GWAS ranged from 0.01% to 6.18%. To further evaluate the potential roles of candidate genes, we used a multi-omic network containing five additional data sets, including leaf and wood metabolite GWAS layers and coexpression and comethylation networks. We also performed a functional enrichment analysis on coexpression nearest neighbors for each gene model identified by the wood anatomical and morphological trait GWAS analyses. Genes affecting cell wall composition and transport related genes were enriched in wood anatomy and stomatal density trait networks. Signaling and metabolism related genes were also common in networks for stomatal density. For leaf morphology traits (leaf dry and wet weight) the networks were significantly enriched for GO terms related to photosynthetic processes as well as cellular homeostasis. The identified genes provide further insights into the genetic control of these traits, which are important determinants of the suitability and sustainability of improved genotypes for lignocellulosic biofuel production.

A High-Performance Computing Implementation of Iterative Random Forest for the Creation of Predictive Expression Networks.

As time progresses and technology improves, biological data sets are continuously increasing in size. New methods and new implementations of existing methods are needed to keep pace with this increase. In this paper, we present a high-performance computing (HPC)-capable implementation of Iterative Random Forest (iRF). This new implementation enables the explainable-AI eQTL analysis of SNP sets with over a million SNPs. Using this implementation, we also present a new method, iRF Leave One Out Prediction (iRF-LOOP), for the creation of Predictive Expression Networks on the order of 40,000 genes or more. We compare the new implementation of iRF with the previous R version and analyze its time to completion on two of the world's fastest supercomputers, Summit and Titan. We also show iRF-LOOP's ability to capture biologically significant results when creating Predictive Expression Networks. This new implementation of iRF will enable the analysis of biological data sets at scales that were previously not possible.

Finding New Cell Wall Regulatory Genes in Populus trichocarpa Using Multiple Lines of Evidence.

Understanding the regulatory network controlling cell wall biosynthesis is of great interest in Populus trichocarpa, both because of its status as a model woody perennial and its importance for lignocellulosic products. We searched for genes with putatively unknown roles in regulating cell wall biosynthesis using an extended network-based Lines of Evidence (LOE) pipeline to combine multiple omics data sets in P. trichocarpa, including gene coexpression, gene comethylation, population level pairwise SNP correlations, and two distinct SNP-metabolite Genome Wide Association Study (GWAS) layers. By incorporating validation, ranking, and filtering approaches we produced a list of nine high priority gene candidates for involvement in the regulation of cell wall biosynthesis. We subsequently performed a detailed investigation of candidate gene GROWTH-REGULATING FACTOR 9 (PtGRF9). To investigate the role of PtGRF9 in regulating cell wall biosynthesis, we assessed the genome-wide connections of PtGRF9 and a paralog across data layers with functional enrichment analyses, predictive transcription factor binding site analysis, and an independent comparison to eQTN data. Our findings indicate that PtGRF9 likely affects the cell wall by directly repressing genes involved in cell wall biosynthesis, such as PtCCoAOMT and PtMYB.41, and indirectly by regulating homeobox genes. Furthermore, evidence suggests that PtGRF9 paralogs may act as transcriptional co-regulators that direct the global energy usage of the plant. Using our extended pipeline, we show multiple lines of evidence implicating the involvement of these genes in cell wall regulatory functions and demonstrate the value of this method for prioritizing candidate genes for experimental validation.

Plant Host-Associated Mechanisms for Microbial Selection.

Plants serve as host to numerous microorganisms. The members of these microbial communities interact among each other and with the plant, and there is increasing evidence to suggest that the microbial community may promote plant growth, improve drought tolerance, facilitate pathogen defense and even assist in environmental remediation. Therefore, it is important to better understand the mechanisms that influence the composition and structure of microbial communities, and what role the host may play in the recruitment and control of its microbiome. In particular, there is a growing body of research to suggest that plant defense systems not only provide a layer of protection against pathogens but may also actively manage the composition of the overall microbiome. In this review, we provide an overview of the current research into mechanisms employed by the plant host to select for and control its microbiome. We specifically review recent research that expands upon the role of keystone microbial species, phytohormones, and abiotic stress, and in how they relate to plant driven dynamic microbial structuring.

Development of a high efficiency integration system and promoter library for rapid modification of Pseudomonas putida KT2440.

Pseudomonas putida strains are highly robust bacteria known for their ability to efficiently utilize a variety of carbon sources, including aliphatic and aromatic hydrocarbons. Recently, P. putida has been engineered to valorize the lignin stream of a lignocellulosic biomass pretreatment process. Nonetheless, when compared to platform organisms such as Escherichia coli, the toolkit for engineering P. putida is underdeveloped. Heterologous gene expression in particular is problematic. Plasmid instability and copy number variance provide challenges for replicative plasmids, while use of homologous recombination for insertion of DNA into the chromosome is slow and laborious. Further, most heterologous expression efforts to date typically rely on overexpression of exogenous pathways using a handful of poorly characterized promoters. To improve the P. putida toolkit, we developed a rapid genome integration system using the site-specific recombinase from bacteriophage Bxb1 to enable rapid, high efficiency integration of DNA into the P. putida chromosome. We also developed a library of synthetic promoters with various UP elements, -35 sequences, and -10 sequences, as well as different ribosomal binding sites. We tested these promoters using a fluorescent reporter gene, mNeonGreen, to characterize the strength of each promoter, and identified UP-element-promoter-ribosomal binding sites combinations capable of driving a ~150-fold range of protein expression levels. An additional integrating vector was developed that confers more robust kanamycin resistance when integrated at single copy into the chromosome. This genome integration and reporter systems are extensible for testing other genetic parts, such as examining terminator strength, and will allow rapid integration of heterologous pathways for metabolic engineering.

Hybridization leads to interspecific gene flow in Sarracenia (Sarraceniaceae).

PREMISE OF THE STUDY: The role of hybridization in plant evolution remains a source of intense debate. Potential consequences range from genetic dead-ends to species fusion or hybrid speciation. While much has been learned from model systems such as Populus, Iris, and Helianthus, many questions remain. Consisting of 11 species that are all capable of hybridizing, Sarracenia presents an excellent system in which to study hybridization. • METHODS: Using microsatellites, we examined a single field site consisting of three species: S. leucophylla, S. alata, and S. rubra subsp. wherryi. We determined the level of genetic admixture and compared it with allopatric sites of the same taxa. • KEY RESULTS: In contrast to the well-defined clusters formed when we examined the allopatric sites, the sympatric field site exhibited a wide range of admixture. Additionally, when the relative genetic makeup of "pure" species at the site was compared with the makeup of hybrids, we found that Sarracenia alata contributed disproportionately to the hybrid genomes. • CONCLUSIONS: Our study provides further evidence that hybridization is contributing to interspecific gene flow in the genus and that all species do not contribute equally to hybridization. Implications for conservation are discussed.

Genetic diversity in three endangered pitcher plant species (Sarracenia; Sarraceniaceae) is lower than widespread congeners.

PREMISE OF THE STUDY: Narrow-ranging, rare species often exhibit levels of genetic diversity lower than more common or widespread congeners. These taxa are at increased risk of extinction due to threats associated with natural as well as anthropogenic events. We assessed genetic variation in three federally endangered Sarracenia species. We discuss maintenance of genetic diversity and evolutionary implications of rarity. • METHODS: We analyzed three noncoding chloroplast regions and nine microsatellite loci in populations spanning the geographic ranges of S. oreophila, S. alabamensis, and S. jonesii. The same microsatellite loci were used to examine a single field site of three more widespread species (S. alata, S. leucophylla, and S. rubra subsp. wherryi). • KEY RESULTS: All three endangered species have experienced reductions in population size and numbers. All show considerably less variation than more widespread members of the genus. Sarracenia alabamensis maintains the greatest microsatellite variation but has the fewest remaining populations and may be under the greatest threat. More widespread S. oreophila maintains surprising chloroplast diversity, yet exhibits little microsatellite diversity. Sarracenia jonesii lacks chloroplast diversity, yet maintains greater microsatellite diversity than S. oreophila. • CONCLUSIONS: The three endangered species differ in levels and structure of diversity, yet not in predictable ways, emphasizing that unique demographic and ecological histories, rather than current distribution and population size, best explain present patterns of genetic variation. Maintenance of remaining genetic variation is important, but preventing further habitat loss and degradation is critical.